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We explore the feasibility to study exclusive meson production (EMP) in hard 
muon-proton scattering fi p — * fi p M at the COMPASS experiment. These mea- 
surements constrain the off-forward parton distributions (OFPD's) of the proton, 
which are related to the quark orbital contribution to the proton spin. 



1 Introduction 

The internal structure of the proton has been a subject of high interest for 
many decades. In general, the spin of the proton resides in its quark and gluon 
constituents and in their orbital angular momenta: 

i(Aq„) + i(A5) + (AG) + L 9 + i 3 = i (1) 

Here, ^{Aq v } denotes the spin in the valance quarks, ^(AS) is the total spin 
carried by the sea quarks and antiquarks, (AG) represents the spin residing in 
the gluonic field, L q and L g represent the orbital momenta of the quarks and 
gluons, respectively. A main goal of the COMPASS experiment at CERN is 
the determination of the gluon contribution AG. In COMPASS, AG will be 
studied via the photon-gluon fusion process, taggeieither by charm- ant icharm 
production, or by high-p^ hadron pair production^. 

A factorizatiorL-theorem for deepLfcvirtual Compton scattering (DVCS) 
was provided by Jicl and Radyushkin an. For exclusive meson production by 
longitudinally polarized photons, a factorization theorem was proven recently 
by Collins, Frankfurt, and Strikmancl. The theorem allows one to separate 
the meson production amplitudes into a perturbatively calculable component, 
describing the interaction of the virtual photon with quarks and gluons of 
the target, and a matrix element component which contains all information 
about the long-distance non-perturbative strong interaction dynamics in the 
produced meson and the nucleon target. The latter component can be param- 
eterized, at leading order PQCD, in terms of 4 generalized structure functions, 
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Figure 1: Illustration of the fac- 
torization theorern-ifor exclusive 
meson production Q. 



Figure 2: Differential cross section for 7r° produc- 
tion at t = t m i n , Q 2 = 4 GeV 2 , as a function of 
xg, in units of nb/GeV 2 , for i th* j l hiwcj different 
valence quark paramctrizationsEijEjliZltj- 



the so called skewed or off-forward parton distributions (OFPD's). Ji showed 
□ that the second moment of these OFPD's determines the total quark angu- 
lar momentum, the sum of the intrinsic quark spins and their orbital angular 
momentum. The OFPD's differ from the usual parton distributions probed 
in inclusive reactions by having a non-zero momentum transfer between the 
proton in the initial and final state. Therefore, the OFPD's are not acces- 
sible in standard inclusive measurements. They can, however, be measured 
in deeply-virtual Compton scattering or in hard exclusive lepto-production of 
mesons (EMP). Thus, exclusive meson production provides a unique tool to 
explore the contribution of orbital angular momentum and the valance quark 
polarization to the nucleon spin. 

EMP are important not only because they relate to the quark orbital mo- 
mentum contribution to the proton spin. They should test our understanding 
of QCD generally, and of OFPD's in particular. COMPASS may lead in the 
study of these new processes, that are only now calculable in QCD. 

In principle, Deeply Virtual Compton Scattering represents the most el- 
ementary process to access the OFPD's. This reaction has been considered 
for COMPASS It was, however, realized that this channel may suffer from 
potentially large backgrounds. Detailed understanding of (exclusive) meson 
production is therefore a precondition for studying DVCS. In turn, DVCS may 
represent an important complement for the Exclusive Meson Production study, 
which we discuss below. 

With the high muon beam energy and high luminosity in COMPASS 
experiment, large Q 2 values can be reached with reasonable rates. For po for 
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example, a Q 2 range between 2 and 20 (80) GeV 2 is accessible for incident 100 
GeV/c (200 GeV/c) muon beam. Thus, this COMPASS program extends the 
approved activity of the CLAS collaboration at JLABtJ, which will explore the 
OFPD in a limited Q 2 range up to about 4 GeV 2 ^We distinguish between EMP 
channels in which the final state meson (p°, 7r + EI, 0, etc.) can be detected via 
magnetic spectrometers, and those which require electromagnetic calorimetry. 
The charged cases can be detected with higher resolution and more easily 
without EM calorimetry. 

As another important aspect of this study, we also consider EMP of a "di- 
pion" (both resonant ^nd non-resonant), such as 7r + 7r~, KK, ttK, as studied 
recently by PolyakovEJ. These works demonstrated that the shape (not ab- 
solute cross sections!) of these "di-pion" mass distributions carry important 
information about quark distribution in the pion (kaon), and also about the 
structure of resonances (p,f%,K*, etc.). 

2 The Proposed Experiment 

2.1 Experimental methodology and trigger 

We base our considerations on the scattered muon trigger that is forseen in 
the COMPASS experiment to study spin dependent DIS. The trigger fires if 
the muon is scattered at large angle 6^ . In our simulations we assumed a good 
angular acceptance for the range 10(5) < # M < 50 mrad for 100 (200) GeV 
beam energy. The trigger is based on the X-Y correlations between hodoscope 
elements at about 31 and 50 meters downstream of the target. 

Average recoil proton energies for EMP are quite low (~30-40 MeV). Thus, 
the recoil proton kinematics for COMPASS do not allow installing an efficient 
proton recoil detector to isolate "elastic" reactions. In order to be able to make 
optimal use of the experimental setup employed by the main muon program, 
or even to run initially in parallel with the main program, we will evaluate the 
inelasticity / oc {M 2 — M 2 ) using measured momenta of the incident and final 
muons and produced mesonEH. The resolution will, however, not be sufficient 
to cleanly resolve events with a recoil proton or a recoil A + or— N* in the 
final state. However, the OFPD's for these processes are related EH. Thus, a 
comparison of the data should be made to model calculations which combine 
both, proton, A + and N* production. 

The main goal of the proposed program is to obtain high statistics EMP 
data on a pure unpolarized proton target during dedicated data runs. This 
would also help to avoid "dilution" by the presence of different target nuclei 
in the polarized proton target. Also, with the liquid hydrogen target, the 
detection efficiency of produced tt° decay 7's would increase due to the bigger 
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total radiation length of the LH. However, in order to explore the feasibility of 
this program already during the initial data taking period of COMPASS, we 
assumed the standard polarized target for our present simulations and count 
rate estimates. 



2.2 Gross section for exclusive meson production (EMP) 

We estimate the dependence of the cross section for the reaction /i p — > /i p ir° 
on the kinematic variables JD 2 , xg, and t, using models for off- forward parton 
distributions. Following Era, we constructed the relevant off- forward parton 
distributions starting from polarized valence quark distributions. 

The exact QCD calculation involves time-consuming numerical integra- 
tion, which is inconvenient for Monte Carlo purposes. However, neglecting 
QCD logarithmic evolution of off-forward parton distributions, one can obtain 
a simple parametrization based on a fit to the results of the exact formulae: 



J(7LP--° p) = ^(x,Q 2 ) 



x e B (*-W0, (2) 



where: 



as 2 (Q 2 )-PF(a:,Q 2 )-UF(x) 



Q 2 (Q 2 +AP) 2 

Here, as(Q 2 ) is the one-loop QCD coupling constant, which - as usual for LO 
calculations - was evaluated with Aqcd = 0.2 GeV and with NF = 4 for the 
number of active flavors. In this equation, VF(x,Q 2 ) is given by: 

PF(x, Q 2 ) = Q4 (1 ~ X) j. (4) 

(-Q 2 - 0.881721 x+ Q 2 x) 2 

UF(x) is the part of the cross-section which is independent of Q 2 . It would 
acquire a logarithmic Q 2 dependence when QCD evolution effects are included. 
The shape of UF(a;) depends on the model of the off-forward parton distribu- 
tions used in the actual calculation. 

Figure 2 shows the differential cross section for 7r° production at t = t m j n , 
Q 2 = 4 GeV 2 , as a function of Bjorkcn x^, in units of nb/GeV 2 , for the 
three different valence quark parametrizations. Thc_solid line refers to the 
Gehrmann-Stirling polarized valence para m fi trization l 15 l the short-dashed line 
is based on the Carlitz-Kaur model E^iiMZl and the loiig-dashed line results 
from the Bartelski-Tatur polarized valence distribution E3. We do not include 
twist- 2 contributions to the OFPD, which according to recent calculations EJEil 
may change the absolute cross sections for EMP. 
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2.3 Monte Carlo results 



The events of exclusive 7r° production, /i+p — > n + p + ir with the subsequent 
decay tt° —> 77, were generated with rates proportional to the estimates of 
differential cross sections from Ref. 18. The directly measured cross section 
could be expressed in terms of the virtual photoproduction cross section: 

d<Vy-> M JV7r°(^Q 2 ,*) _ rfn 2 ^ doyjv^jy^oO = Q 2 /(2M p v),Q 2 ,t) 

d^dt " 1 [Q ' V] dt ' (5) 

where T is the flux of virtual photons: 



(6) 



Here, a is the fine structure constant, and e is the virtual photon polarization. 

The cross sections were estimated for proton and carbon targets. The latter 
one is more representative for the complex material of the COMPASS polarized 
target. For the proton target, der 7 * n^n-k° (x, Q 2 , t) /dt is given by the Eqs. (||- 
[|) and B — 4.5 GeV~ 2 . For x = 0.1, Q 2 = 4 GeV 2 and t = t ml — % is equal 
to 0.88 nb/GeV 2 . For similar conditions, Vanderhaeghen et al. 113 give about 
4 times larger cross section. We will therefore give counting rate estimates 
in this report with an uncertainty, which reflects the difference between these 
two theoretical estimates. The present theoretical uncertainties are mainly due 
to different assumptions on the shapes/magnitudes of the valence quark spin 
distributions. For the carbon target, the cross section was assumed to be the 
product of the proton cross section and a function: 

/(*) = 8.12 e - 4 7.7(t-W) + 0.3 . (7) 

This function was obtained from the parametrization of the NMC data for 
exclusive p° production on a carbon targetOJl 

In our simulations, we considered the kinematic range of 1 GeV 2 < Q 2 < 
Q 2 nax GeV 2 , and 10(15)GeV < v < 85(180) GeV. Here, Q 2 max = AE^E^ - v) ■ 
sin 2 (0™ aa: /2), Q" lax — 50mrad , and is the incident muon energy. The values 
outside (inside) the brackets correspond to the beam energy of 100 (200) GeV. 
The same convention will be used hereafter. The upper cut on v was chosen to 
eliminate the kinematic region where the amount of radiative events is large, 
whereas the lower one to eliminate events with poor kinematic resolution. The 
value of Qmax is about 20 (90) GeV 2 . The cross section integrated within this 
kinematic regime 
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Figure 3: The total acceptance to observe exclusive n° events as a function of v (top) and 
Q 2 (bottom) at a beam energy of 100 GeV. The different histograms correspond to different 
geometrical acceptances of the electro-magnetic calorimeter 

The total acceptance e^o to observe an exclusive ir° event, taking into 
account the geometric acceptance, the absorption of the decay photons before 
entering the calorimeter, and the reconstruction of tt° from the calorimetric 
measurements, is shown in the top parts of Fig. 3 as function of v by the 
solid histograms. The main effect on e^o is due to a poor reconstruction of low 
energy ir° mesons from the calorimetric measurements, at small v, and due to 
the absorption of photons in the polarized target, most significantly at large 
v. The average total efficiency e^o is 0.43 for both incident beam energies. For 
the kinematic range v > 15 GeV, Q 2 > 1 GeV 2 , 0.02 < x < 0.2, the efficiency 
is quite large, varying between 0.3 and 0.5. 

Taking the above values of cross sections, and the expected luminosity 
for the COMPASS experiment!] of 43 pb _1 day _1 , we evaluate the number of 
deep inelastic exclusive 7r° events per day, N^ 0N /day(Q 2 > 1 GeV 2 ), to be 
about 248-992 (53-212) for 100 (200) GeV. After taking into account effects 
of acceptance, secondary interactions, and the efficiency of tt° reconstruction, 
the number of events, Npf N /da,y(Q 2 > 1 GeV 2 ), is about 107-428 (23-92). 
The expected numbers of exclusive ir° events for a running period of 150 days 
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(corresponding to a 1 year period) assuming an overall (SPS and COMPASS) 
efficiency of 25% are 4013-16050 (855-3420). For Q 2 > 4GeV 2 , we expect only 
130 (42) counts. Adopting the Carlitz-Kauer valence quark parametrization 
the last numbers will increase by about a factor of 4. Thus, in 1 year of 
COMPASS operation we will collect at most 500 tt° events with Q 2 > 4GeV 2 . 

Up to this point we assumed full size electro-magnetic calorimeters (total 
of 18 m 2 ) at their default positions. The different curves in Fig. 3 show the 
result if only limited coverage is assumed. While a reasonable efficiency may 
still be obtained with about half the calorimeter size (9 m 2 ), a further cut (e.g. 
down to 7 m 2 ) significantly reduces the ir° efficiency. Furthermore, the count 
rates presented above were based on a polarized target. Assuming that the 
overall length of a target is limited to typically 2 m, a dedicated LH target 
would lead to a reduced luminosity (about a factor of 4) unless the beam 
intensity can be correspondingly increased. 

The detailed discussion of the possibility for COMPASS experiment to 
measure another exclusive vector meson production channel, fi + N — > /i + 
N + p° , is presented in Ref. c3. Here we only briefly summarize the results 
of the Monte Carlo calculations for this reaction. The kinematic range is 
somewhat narrower than for the previous channel: 2GeV 2 < Q 2 < Q 2 max and 
20(40) < v < 90(180) GeV. The estimated value of the measured cross section 
is equal to 0.46 nb (0.5 nb) for 100 GeV (200 GeV). These values are 
obtained using a parametrization of the cross section measured in the NMC 
experiment cB. For the given kinematic range, the average total acceptance 
e p o is 0.73 (0.68). Taking into account the expected luminosity for the COM- 
PASS experiment, and the effects of acceptance and secondary interactions, 
the number of events with the invariant mass of two pions in the central part 
of p" peak, 0.62 < m vv < 0.92 GeV,A^yday(Q 2 > 2 GeV 2 ), is about 10 4 
for both beam momenta. For a running period of 150 days assuming again an 
overall (SPS and COMPASS) efficiency of 25% this will result in about 3.7-10 5 
events. For Q 2 > 4GeV 2 we expect 75 k (107 k) events. 

3 Conclusions 

We explored the feasibility to study exclusive meson production in hard muon- 
proton scattering p — ► fj, p M at the COMPASS experiment. For moderate 
Q 2 > 1 GeV 2 , a significant counting rate for exclusive 7r° production can be 
reached. For the experimental conditions assumed in our calculations, still 
up to 500 events beyond a Q 2 of 4 GeV 2 can be collected during 1 year of 
operation. 

On the other hand, COMPASS may provide already during the first main 
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H runs, high statistics data on exclusive p° production even at Q 2 beyond val- 
ues presently accessible at JLAB. Furthermore, data taken with the polarized 
target may also allow a first view at polarization phenomena. 
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